Phase Drifts of Sub-pulses during the 2004 Giant Flare of SGR 1806-20 and 

Settling of the Magnetic Fields 



Yi XING and Wenfei YU 

Key Laboratory for Research in Galaxies and Cosmology and Shanghai Astronomical Observatory, 
80 Nandan Road, Shanghai 200030, China. E-mail: wenfei@shao.ac.cn 

ABSTRACT 

We analyzed the observations of SGR 1806-20 performed with the Rossi X-ray Tim- 
ing Explorer (RXTE) during its 2004 giant flare. We studied the phase evolution of 
the sub-pulses identified in the X-ray waveform and found that the sub-pulses varied in 
phase with time and then gradually settled, which might indicate drifts of the emission 
regions in relative to the neutron star surface, or changes in the local emission geom- 
etry before the magnetic field became stable. The characteristic e-folding timescale of 
the phase drifts measured starting about 15 s following the initial flux spike are in the 
range between 37 s and 84 s. This leads to the first measurements of the characteristic 
timescale for the magnetic field of the neutron star to settle after a field reconfiguration 
during the giant flare. 

Subject headings: neutron stars, Soft Gamma-ray Repeaters 



Introduction 



Soft Gamma-Repeaters (SGRs) are thought to be a type of magnetars: isolated neutron stars 
with much higher magnetic field than normal pulsars. The characteristic strength of magnetic 
field is believed to be about 10^^ — 10^^ G, inferred from the observational measurements of the 
pulse period and its derivative assuming a dipole field. Recurrent X-ray and soft gamma-ray flares 
have been seen from SGRs, with a typical luminosity up to lO^'^ ergs/s since their discovery in 
1979 (for an observational review, see Mereghetti 2008). Several major observational properties of 
the flares in SGRs can be explained by the most popular model in which the flares are produced 
by Alfven wave which propagates in the magnetosphere of the neutr on star, due to the magneti c 
field diffusion associated with small cracks in the neutron star crust (iThompson &: Duncanlll995l ). 
Besides these recurrent flares, several giant flares with much higher luminosity have been observed. 
The luminosity during giant flares could reach more than a million times the Eddington luminosity 
of a common neutron star, i.e. 10^^ — 10^^ erg/s. As an example, for the most extre me giant flare 



observed from SGR 1806-20 in 2004, the peak luminosity was about 2 x 10^^ erg/s (jHurlev et al 



20051 : iMereghetti et al.ll2005l : iPalmer et al.ll2005l : iTerasawa et al.ll2005l ). Energy of giant flares was 
thought to have a magnetospheric origin other than a neutron star crust origin. In the above 
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magnetar model, giant flares are caused by a fireball of pair-dominated plasma that expands at 
relativistic speed. The unwinding internal magnetic field will build up the magnetic energy outside 
the neutron star slowly by pushing out inside electric currents into the magnetosphere. When the 
magnetosphere currents reach a point at which the magnetosphere becomes dynamical unstable, 
the twist field outside the star associated with dissipation will have a sudden relaxation and the 
magnetic field will have a reconfiguration, consequently the destructive magnetic instability may 
produce the neutron star crust fracturing in large scale, which can lead to an extreme outburst 
(jxhompson k DuncanlflQQsl : iLvutikov l[2006l l. 

The total energy inferred from the magnetic strength of magnetar is of about 10^^ erg, which is 
signific antly more than energy conta ined in common neutron stars in the form of rotational kinetic 
energy (jBraithwaite &: Spruit II2006I ). During a giant flare a large portion of the energy contained 
in magnetar would be lost, so the giant flares are very rarely seen. Up to now only three events have 
been detected: the 1979 flare from SGR 0526-66, the 1998 flare from SGR 1900+14, and the 2004 
flare from SGR 1806-20. Interestingly, it was suggested before the 2004 Giant flare that magnetars 
would have a supergiant flare produced by a dynamic overturning instability. These type of flare 
should be 100 times less frequent than normal giant flares bu t could be obs erved within about 10 
times the distance because of its much higher energy release (jEichler II2002I ). The 2004 giant flare 
was just like the predicted events. 

Giant flares have some common properties. In the X-ray waveband the light curves of giant 
flares consist of a short hard initial spike followed by a long softe r pulsating tail. The peak luminosi- 



ties of the 1979 flare and the 1998 flare were abo ut 10^^ erg/s (IMazets et al.lll979l : iFenimore et al 



1981 



Hurlev et al 



1999; 



erg/s (jHurlev et al.l l2005l : iMereghetti et al. 



Feroci et al 



19991. 200lh while during the 20 4 flare it reached to 2 x lO-^^ 



2005 



Palmer et al. 



2005 



Terasawa et al.l 120051 ). The 



rise time is typically smaller than a few milliseconds during which the luminosities of giant flares 
increase rapidly. After the initial spike there is the pulsating tail lasting about hundreds of seconds, 
which is thought to originate from the cooling of the ejected plasma which remains conflned in a 
trapped thermal flreball near the stellar surface due to the strong magnetic fleld in magnetar model 
(jThompson &: Duncanlll995l ). 



For the three known giant flares, the pulsating tails always show a strong periodic modulation 
at the neutron star spin period and have complex pulse proflles. They are usually not single- 
peaked and show strong evolution with time. During the 1998 flare from SGR 1900+14, the pulse 
profile of the pulsating tail showed four p eaks and a dip wh ich were almost evenly spaced at about 

and the four-peak structure evolved 



1.03 s interval in each period of 5.16 s (iFeroci et al 



2001 



to being nearly sinusoidal with time. The multi-peak features are thought to be described as 
collimated X-ray jets from trapped fireballs, in comparison to the 1979 giant fiare from SGR 0526- 
66 during which only two sub-pulses have been observed. It was suggested that the number of 
sub- pulses detected in each perio d of these two sources could be due to different number of X-ray 
jets ( Thompson fc: Duncan IlioOll ). For the 1998 giant fiare the phases of these four sub-pulse peaks 
were recognized stable, which may suggest that the X-ray jets are tied to surface features on the 
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neutron star. The four-peaked pa ttern is thought to be direct evidence of multipolar structure of 

Investigations in the quiescent emission of 



the neutron star's magnetic field (jFeroci et al 



2001 



SGR 1900+14 revealed th at before the 1998 g ;iant flare the pulse profile of the source always showed 
complex multi-structure (jWoods et alj|200ll ). and once it evolved to sinusoidal during the flare it 



keeped unchanged for more than 1.5 years. It is suggested that if the multipeak profile indicates 
multipolar field structure and the simple sinusoidal profile indicates dipolar field geometry, during 
the time of the 1998 giant flare the magnetic field, both inside and outside the neutron star, 
underw 'ent a global reconfi guration, while during the quiescent stage the magnetic field was relative 



stable (Woods et al. 2001 



The evolution of pulse profile was also observed during the 2004 flare from the SGR 1806- 
20. It showed an opposite trend compared with the 19 98 flare; the pulse profile e volved from a 
simple sinusoidal to a complex multi-peaked structure (jWoods Sz Thompson! l2006l ) . The phases 
of these multi-peak ar e generally though t stable, indicatiiig that there was no subs tantially field 



change during the tail ([Palmer et al.ll2005l : iMereghetti et al. 



2005 



Boggs. et al. 1120071 ) , even though 



we see obvious phase shifts of the sub-pulses on short time scales (see our study below). In the 
quiescent X-ray light curve prelude to the 2004 flare, the pulse profile was nearly simple sinusoidal, 
it became more complex after the flare. Considering the pulse proflle changes during the flare, 
it was recognized that the pul se profile changes in the quiescent light curve were fiare-induced, 
consistent with SGR 1900+14 (jWoods et aLlboOTt ). 



Besides the two SGRs from which evolving pulse profile has been observed, during magnetar 
flares from magnetar candidates, for example CXOU J164710. 2-455216, the pulse proflle was also 
observed to evolve. During the observations of 4.3 day prior to and 1.5 day subsequent to two 
remarkable events that were detected with Swift on 2006 September 21 (a 20-ms burst and a 



glitch ) , the pulse proflle was found to change from single-peaked to three distinct peaks (jMuno et al 



20071 ). For this type of source with lower X-ray energy release, a large-scale rearrangement of the 
magnetic fleld was thought unreasonable. It is suggested that the pulse shape was governed by 
the distribution of currents within the magnetosphere, and t he proflle change s hould be due to the 
change in the distribution of currents in the magnetosphere (iMuno et al.l 120071 ). 



In the following study we focus on the pulse proflle during the pulsating tail of the 2004 flare 
form SGR 1806-20. In our previous study of SGR 1806-20 giant flare using GoodXenon mode 
data of the RXTE/PCA only, we have decomposed this pulse proflle into four sub-pulses. We 
modeled the pulse proflle with four Gaussian functions and a constant background flux level, and 
investigated the X-ray waveform evolution o f each sub-pulses, excluding the flrst 2-3 cycles because 
of data gaps in the GoodXenon mode data (jXing &: Yull2009l ). In that study we found the phases 
of the sub-pulses were not flxed during the tail. The phase of one sub-pulse was decreasing while 
the phase of another sub-pulse was increasing with time - an apparent anti-correlation between the 
phase drifts of the two sub-pulses. This suggests that the phases of the sub-pulses moved although 
the multi-peak waveform structure remained unchanged. Modeling the sub-pulses as multi-peak 
Gaussian gives good description of the overall phase distribution of the photons in each sub-pulses. 
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but not accurate and sensitive in determination of the phases of the sub-pulse peaks. Furthermore, 
the data gap in the GoodXenon mode data recorded at the beginning of the Giant flare because of 
telemetry saturation prevents us from studying the sub-pulse phase evolution during the first few 
cycles of the giant fiare. In this study, we intend to measure the phases of the sub-pulse flux peaks, 
which show the sub-pulse phase shifts more obviously. We determined the sub-pulse phase by 7- 
point Gaussian fit to each sub-pulse peak with Standard 1 mode data instead of GoodXenon data. 
This data mode does not have data gaps and uniquely helps determination of the sub-pulse phases 
during the first few cycles. We found strong evidence of systematic phase shifts of those sub-pulse 
peaks, in agreement with our previous results but further helped us measure the characteristic 
timescales of the phase evolution of the sub-pulses. This timescale might indicate the timescale for 
the magnetic field to settle on the neutron star after a reconfiguration, which has been suggested 
to associate with the giant energy release during the flare. 



2. Observations and data analysis 

The giant flare of SGR 1806-20 was recorded by the Rossi X-ray Timing Explorer (RXTE) on 
2004 December 27th, during that time the RXTE was pointing at the cluster of galaxies A2163 (with 
the RXTE Observation ID: 90132-01-13-09). The data we used was taken by Proportional Counter 
Array (PGA) in GoodXenon mode with a time resolution of 9.5367431640625e-07 s (for determining 
pulse period and checking if the sub-pulse phase drift was not obviously energy dependent) and 
Standard-1 mode data with a time resolution of 0.125 s (for investigating sub-pulse phase drift). 

We used the Standard-1 mode data to generate the light curve of SGR 1806-20 during the 
giant flare in Figure 1. The time resolution is 0.125 s. The X-ray emission during the giant flare 
consists of two components: a descending non-pulsed component and a pulsed component. In order 
to accurately determine pulse period and sub-pulse phase during the tail, we need to remove the 



non-p ulsed component from the light curve. This was not done in our previous work (IXing &: Yu 



20091 ) . We separated the light curve into segments of cycles starting at the flrst period during the 
giant flare (the phase-zero was selected at the start of the 3rd period, corresponding to the RXTE 
raw space craft c lock time of 346800645.87843037 s) using a rough period of 7.56 s as reported 



([Israel et al.ll2005l ) and calculated the average count rates for each cycle. This is shown in Figure 
1. 3-points quadratic interpolation was applied to the average count rates in order to construct 
the non-pulsed component. By removing the non-pulsed component from the high time resolution 
light curve (1/4096 s, generated from GoodXenon mode data), we obtained the pulsed component 
during the pulsating tail, which was used to determine the pulsed component and the sub-pulse 
phase shift. Gonsidering the rough period was determined in a way influenced by signiflcant proflle 
changes during the first few cycles, we determined the pulse period during a 400s interval starting 
at about 49 seconds after the initial spike when the pulse profile was relatively stable. We folded the 
400 s data into 128 phase bins and searched for the best period in the period range from 7.5536 s to 
7.5663 s using a step-length of 0.0001 s. A best period of 7.5611 s was obtained, which is different 
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from the 7.5605 s period determined in our previous approach without removing the non-pulsed 
component. Using this period, we generated pulse waveform of each cycle during the pulsating 
tail. We also checked the waveforms in the two energy ranges in which the total counts were nearly 
equal, i.e. 2-15 keV and 15-60 kcV. We concluded that the waveform during the pulsating tail did 
not show obvious energy-dependent phase shift. Therefore we went further using the light curve in 
the entire PCA band to study the sub-pulse evolution. 

We then apply the newly determined pulse period to the Standard-1 mode data to study the 
waveform evolution during the pulsating tail, because the data covers the gaps in the GoodXenon 
mode data at the beginning of the giant flare introduced by the saturation of RXTE telemetry while 
the time resolution is good enough for sampling on the time scale of the pulse period. Following 
the previous procedure to remove the non-pulsed component, we obtained the pulsed component 
during the pulsating tail using the period of 7.5611 s. Then we plotted the waveforms of the 
pulsed component of each cycle in Figure 2. A total of 35 cycles were shown. A portion of the 
inferred pulsed emission during the first cycle was invisible (dashed line) because the non-pulsed 
component was overestimated assuming a linear interpolation. From Figure 2 we could see that 
the overestimation was introduced by the influence of the initial spike in the flrst cycle. 

Multi-Peak structures could be clearly seen from these pulse profiles, as shown in Figure 2. 
Each pulse profile consists of four sub-pulses, and the phases of the sub-pulses were not fixed during 
the fiare. A circle was used to mark the peak value of the major sub-pulse of the first cycle, which 
seemed to relate to the 3rd sub-pulse. This indicates that the waveform evolved drastically at the 
beginning of the giant flare. In order to determine the phases of the sub-pulses, we selected seven 
data points centered at the peak of each sub-pulse and fltted them with a single Gaussian function. 
Notice that these pulse profiles were obtained by removing the trend of the non-pulsed component 
and the mean count rate level during the cycle has been included. We found in most cases the 
sub-pulse peaks could be fitted by Gaussian functions with reduced iii the range between 0.20 
and 9.4. As an example, the pulsed profile of the 8th cycle and the Gaussian components resulted 
from those seven-points Gaussian fits are shown in Figure 3. 

3. Results 

3.1. Phase drifts of sub-pulses 

We show the evolution of the phases of the sub-pulses in Figure 4. We called the sub-pulse 
at about phase 0.2 as the 1st sub-pulse, the sub-pulse at about phase 0.45 as the 2nd sub-pulse, 
the sub-pulse at about 0.55-0.6 as the 3rd sub-pulse, and the sub-pulse at about 0.73 as the 4th 
sub-pulse. The 3rd sub-pulse was excluded in our analysis because it overlaps with the 2nd sub- 
pulse and its phase can not be well-determined. In addition the reduced chi-square values were 
greater than 10 in the Gaussian fits for a few cycles. The large reduced chi-squares indicated that 
the peaks of these sub-pulses can not be modeled by single Gaussian functions. Therefore we did 
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not measure the phases of these sub-pulses, respectively. These cycles are: the first two cycles, the 
16th, the 17th, and the 34th cycles for the 1st sub-pulse; the first two cycles for the 2nd sub-pulse; 
and the 1st, the 16th, the 20th, and the 31th cycles for the 4th sub-pulse. 

Figure 4 shows clearly that the phases of the sub-pulses were not fixed. The evolution trend of 
the 1st sub-pulse and the 4th sub-pulse is opposite, especially in the early stage where the phase of 
the first sub-pulse decreases while the the phase of the fourth sub-pulse increases. We calculated 
the correlation coefficient between these two sub-pulses and got - 0.72ifc0.12, confirm ing that the two 
sub-pulse phase drifts are opposite. As suggested previously in (jXing k. Yull2009l ). the correlation 
may suggest the phase drifts are physically connected. For example, the two sub-pulses may be 
associated with opposite magnetic poles since their phase separation is around 0.5. 



3.2. Characteristic timescales of phase drifts 

We have measured the phase evolution of the 1st sub-pulse since the third cycle and the 4th 
sub-pulse since the second cycle with exponential models to determine the e-folding timescales for 
the phases to become stable. The e-folding time scales are 36.95ib3.01 s for the 1st sub-pulse and 
65.41ib5.99 s for the 4th sub-pulse (see Figure 5). In order to investigate whether smaller e-folding 
timescales would be obtained if we include cycles at earlier stages, which would suggest that our 
results be biased because of lack of data during the first two cycles, we investigated if smaller 
e-folding times would be obtained if we make use of data from earlier cycles. By including cycles 
from the 10th and backward one by one, we checked if the e-folding times obtained from the fits 
would get smaller. We found in general these e-folding times are within uncertainties. Therefore 
the characteristic time scales of tens of seconds we obtained indeed represents the time scale for 
the phases of the sub-pulses to become stable. 

The phase separation of the 2nd sub-pulse and the 3rd sub-pulse seemed more obviously 
changing in the early stage too (see Figure 2), while in the later stage their phases seemed relative 
fixed. The phase of the 2nd sub-pulse showed no gradual change while the phase of the 3rd sub- 
pulse modeled with Gaussian compon ents in our previo us preliminary study was increasing with 



time. Adopting the method used in (IXing &: Yul l2009l ) by decomposing the waveform into four 



Gaussian components, but without a constant component since the long-term background trend 
was removed, we obtained the phase evolution of the 3rd sub-pulse (note: the method is different 
from those for the 1st and the 4th sub-pulses). We determined the e-folding timescale as 84.17±7.54 
s, on which the phase became stable (also see Figure 5). The three best-fit exponential models were 
over-plotted as thick lines in Figure 2 . 

From Figure 2 we could see sub-pulse 1 and 4 are sitting on either sides of sub-pulse 2. Flux 
contribution from sub-pulse 2 may affects our estimates of the phases of sub-pulses 1 and 4, which 
tends to give sub-pulse phases toward sub-pulse 2. In order to investigate this, we also included 
data points of the entire sub-pulses rather than only 7 data points of sub-pulses 1 and 4 as well 
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as the data of the slope on either side of sub-pulse 2. For sub-pulse 4 the minimum phase chosen 
corresponds to the flux minimum during the rising side of the sub-pulse, which was in the phase 
range between 0.5 and 0.7, and the maximum phase chosen was phase 1.0. For sub-pulse 1 the 7 data 
points we selected before and a background data slope covers the entire sub-pulse, i.e., we included 
data points between phase 0.0 and maximum phase of the 7 data points. We then fitted the data 
associated with sub-pulses 1 and 4 to a Gaussian function plus a linear component, representing the 
sub-pulse peak and the sloping side of sub-pulse 2, respectively. The slope of the linear component 
was set to be above for sub-pulse 1 and to be below for sub-pulse 4, as sub-pulse 1 sits on the 
left shoulder of sub-pulse 2, and sub-pulse 4 sits on the right shoulder. An example of the fits for 
the 3rd cycle of sub-pulse 1 and 4 is shown in Figure 6 and 7. We then obtained the phase evolution 
of the sub-pulses. Similarily we ignore those fits with reduced-x^ larger than 10. We then fitted 
the obtained phase evolution with the exponential model to measure the e-folding timescales for 
the phases of the sub-pulses to settle and got 34.05 -|-/- 3.07 s for sub-pulse 1 and 50.94 +/- 6.98 s 
for sub-pulse 4 (see Figure 8). This indicates that the e-folding time scales are indeed on the time 
scales of tens of seconds. 



3.3. Comparison to the decay timescale of the giant flare 



The e-folding decay time of the 1998 giant flare from SGR 1900-1-14 was found to vary with 
energy, e.g., 78 s in the Ulysses ene rgy range of 40-100 keV and 70 s in the Bep poSAX energy range 
of 100-700 keV dFeroci et ahlboOlh . An even longer decay timescale of 90 s (jMazets et al.l 1 19991 ) 
was obtained in the Konus energy range. It was theref ore inferred that t he e-folding decay time of 
the 1998 flare decrease with increasing photon energy (jFeroci et al.ll200ll ). For the 2004 giant flare 



from SGR 1806-20, an e-folding decay timescale of 128.21±0.19 was obtained in the RXTE/PCA 
energy range (2-60 keV). Therefore, the e-folding decay time of the X-ray flux of the giant flare is 
a little larger than the characteristic timescale for the phases of the sub-pulses to settle, but on the 
same order of magnitude. 



4. Conclusion and discussion 



We have determined the phases of the sub-pulses in the X-ray waveform during the 2004 giant 
flare and studied their evolutions. There existed phase drifts of these sub-pulses up to 0.06 in 
less than 100 seconds. The phases of the sub-pulses became stable after a characteristic timescale 
of tens of seconds. This indicates that the magnetic instability responsible for the magnetic field 
reconfiguration damps out on the above timescale and the magnetic field became stable after such 
a characteristic time scale. Two sub-pulses separated by about 0.5 phase, namely the 1st sub-pulse 
and the 4th sub-pulse in our study, show opposite evolution trend. The apparently correlated 
evolution trend is more obvious during the early stage of the giant flare. If the two sub-pulses 
correspond to the emission regions of two opposite magnetic poles or of the same magnetic pole, 



- 8 - 



the phase evolution and the opposite phase trend indicate the two emission regions are probably 
tied together during the giant flare. The ~ 0.5 phase separation of the two sub-pulses and the 
potentially physically correlated phase drifts then likely provide evidence that we can see at least 
two poles of the neutron star. This provides an additional evidence for a multipolar magnetic field 
geometry other than that discussed in Feroci et al (2001) on the four sub-pulse profile. 

The initial spike of a giant flare has a sharp rise which is typically smaller than a few millisec- 
onds. For the 1979 event the rise time was smaller than 2 ms. For the 1998 event it was smaller 



than 4 ms (jMazets et al.lll999l ). During the 2004 giant flare fr om SGR 1806-20 th e rise time was 
much smaller than ever observed, which took only about 0.3 ms (jPalmer et al.ll2005l ). The steep rise 
time of giant flare is thought associated with the Alfven crossing time of the inner magnetosphere, 
which is as short as 30//s estimated by RnsIc, where Rns is the radius of the neutron star and c 
is the light speed (ILyutikov II2006I : iTanaka et al.ll2007l ). During the steep rise the twist field outside 
the neutron star suffers a sudden relaxation, which leads to a rapid release of energy stored in the 



magnetosphere (jLyutikov 



2006 



Braithwaite &: Spruit 1120061 ). Besides the initial steep rise there 



was another timescale noted in previous studies, i.e. the intermediate rise time to the peak. The 
intermediate rise time to the peak during the giant fiare from SGR 1900+14 was 3.1 ms, which 
was shorter than the 9.4 ms observed during the 2004 flare from SGR 1806-20. This timescale is 
thought to be limited b y the propagation o f a fracture, so it could be used to estimate the fracture 
size during giant flares (jTanaka et al.l 120071 ). 



For all of the three giant flares observed, the duration of the initial spike was about hundreds 
of milliseconds: 0.25 s for the 1979 giant fiare from SGR 0526-66; 0.35 s for the 1998 giant flar e 
from SGR 1900+14; and 0.5 s for the 2004 giant flare from SGR 1806-20 (|Mereghetti I bood i . 
It may reflect the Alfven wave crossing time of the star which is in the range between 0.2 s 
and 0.5 s estimated by Rms /Va.ns, where Va,ns is the Alfven wave crossing speed in the star 
(iThompson &: Duncanlll995l ). 



The e-folding decay time of the giant flare was thought to describe the cooling time of the 
trapped flreball, which was formed by e nergy released in the initial s pike trapped in the neutron star 



magn etosphere by high magnetic field (IThompson &: Duncan 



1995 



Palmer et al. 



2005 



Mereghetti 



20081 ) . The emission from trapped flreball in higher energy band seemed to cool a little faster than 
that in lower energy band, which leads to a softening of the dominated emission during the later 
part of the pulsating tail. The characteristic time scales of the phase drifts of the sub-pulses are a 
little smaller, but comparable to the decay time scale of the giant flares. We therefore infer that the 
magnetic field reconfiguration settle earlier than the cooling of the trapped flreball in the magnetar 
model. 

In this work we have found the phases of the sub-pulses were obviously drifting during the 
flrst ten cycles or so and gradually became stable. The characteristic timescale of phase drifts is in 
the range 37 s - 84 s. This may tell us the instability of the magnetic field during the giant flare 
damped in tens of seconds before the held conflguration became stable. Thus this time scale may 
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be used to constrain the magnetic fields and the mechanism for the energy release of the giant flare. 
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Fig. 1.— Top: the (2-60 kcV) RXTE/PCA light curve during the 2004 outburst of SGR 1806- 
20 generated from the Standard- 1 mode data. The time resolution of the light curve is 0.125 s. 
Bottom: evolution of the average rates of each cycles during the pulsating tail. 
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Fig. 2. — Pulse profile evolution during the first 35 cycles of the 2004 giant flare obtained from 
Standard 1 data. The non-pulsed background trend is removed from the data. Time sequence 
increases from bottom to top. Count rate corresponding to each cycle was shifted by 1500 counts/s 
one by one. The circle marks the peak of the major sub-pulse during the first period, consistent 
with the sub-pulse #3 identified in subsequent cycles. 
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Fig. 3. — An example of the model fit to the waveform in which each of the three major sub-pulse 
peaks were fitted with a Gaussian (seven-point fits). The best-fit Gaussians were over-plotted as 
dashed lines. The central phases of the Gaussians gives the phases of the sub-pulses. A constant 
level at the mean count rate of the cycle has been added. 
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Fig. 4. — Phase evolution of the sub-pulses #1, #2, and #4 measured by seven-points Gaussian 
fits to the sub-pulse peaks. 
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Fig. 5. — Phase evolution of the sub-pulses modeled with exponential phase drift. Notice that the 
phases of sub-pulse #3 was obtained from a decomposition of the entire waveform into 4 Gaussian 
components and the phases of the #1 and #4 were obtained from 7-points Gaussian fits. The 
phase drifts should be considered as those in relative to the sub-pulse #2 - the major sub-pulse. 
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Fig. 6. — The model fit to data around sub-pulse 1 considering the left shoulder of sub-pulse 2 of 
the 3rd period. 
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Fig. 7. — The model fit to data around sub-pulse 4 considering the right shoulder of sub-pulse 2 of 
the 3rd period. 
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Fig. 8. — Phase evolution of the sub-pulses modeled with exponential phase drift. Different from 
Figure 5, the phases of the sub-pulses 1 and 4 were obtained from fitting with the sloping sides 
of sub-pulse 2 to a Gaussian function plus a linear slope. Notice that data with reduced-x^ larger 
than 10 were excluded in the plot. Around period No. 17-23 the model does not fit sub-pulse 4 
well. 
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